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Abstract 

RF  power  performances  of  GaN-MESFETs  are  reported  using  a physics-based  model  that 
incorporates  dispersion  in  the  output  resistance  and  transconductance  due  to  traps  and 
thermal  effects.  Calculated  I-V  characteristics  are  in  excellent  agreement  with  the  measured 
results.  Taking  thermal  effects  into  account  the  maximum  output  power  of  a 0.3  pan  x 100/an 
GaN  MESFET  is  22dBm  at  a power  gain  of  4.2dB  at  4GHz.  The  corresponding  quantities  are 
27dBm  and  6.4dB,  respectively  if  a constant  channel  temperature  of  300K  is  assumed.  At 
elevated  temperatures  compression  in  output  power,  gain  and  PAE  is  less  in  MESFETs  with 
longer  gate  lengths. 


Introduction 

Currently,  power  transistors  used  in  applications  involving  phased  array  radar  and 
wireless  base  stations  are  Si-LDMOS  which  operates  around  lW/mm  power  density  range. 
Recently,  GaN-based  FETs  have  successfully  demonstrated  their  usefulness  in  circuits 
operating  at  high  temperatures,  high  power  and  high  frequencies  [1-4].  GaN  with  bandgap  of 
3.4eV,  saturation  velocity  of  3xl07  cm/s  and  low  parasitics  is  suitable  for  high  power  and 
high  frequency  applications  compared  to  Si-LDMOS  and  SiC-based  transistors.  Moreover, 
GaN  grown  on  SiC  offers  a thermal  conductivity  of  4.5  W/cm/K  making  the  system  suitable 
for  applications  at  high  temperatures  and  high  power.  GaN  HEMTs  with  output  power  of 
11.7  W/mm  at  10GHz  with  43%  power  added  efficiency  [1  ],fr-  107  GHz  [2]  and  fmax  = 155 
GFlz  [3],  and  operating  temperature  of  750C  [4]  have  been  reported.  However,  the  inability 
to  dissipate  the  generated  heat  leads  to  self-heating  effects  [5-7].  For  a GaN  HFET  grown  on 
SiC,  with  gate  width  of  250/an,  Nuttinck  et  al.  [5]  estimated  the  maximum  temperature  to  be 
96C  in  the  active  channel  with  the  device  operating  under  a continuous  wave  dissipation 
power  density  of  5 W/mm.  Due  to  thermal  effects  a 25%  reduction  in  drain  current  was 
observed  as  the  drain  pulse  width  increased  from  1%  to  100%  for  the  same  250 pan  wide  GaN 
HFET  grown  on  SiC  [6].  Using  Raman  spectroscopy  Kuball  et  al.  [7]  estimated  the 
temperature  at  the  gate-drain  opening  to  be  as  high  as  180C  for  a Apan  x 200 pan  GaN  HFET 
grown  on  Sapphire  at  drain  and  gate  bias  of  20V  and  0V,  respectively. 
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Self-heating  effects  at  large  drain  bias  increase  the  device  lattice  temperature  and  reduce 
the  transport  parameters  such  as  mobility  and  carrier  saturation  velocity  by  increasing  the 
carrier-phonon  scattering  [8,9].  Besides,  GaN  based  devices  are  plagued  by  traps  that  results 
in  performance  deteriorating  gate-  and  drain-lag  transients  and  current  collapse  in  the  I-V 
characteristics  [10,  11].  Trapping  effects  are  dependent  upon  the  applied  signal  frequency 
and  their  effect  is  more  pronounced  for  frequencies  less  than  the  transconductance  and  output 
resistance  dispersion  frequencies  [12,13].  The  dispersion  frequencies  are  related  to  the 
detrapping  time  constants,  which  are  critically  dependent  upon  the  junction  temperature. 
Therefore,  an  exact  simulation  of  GaN  based  devices  should  proceed  by  determining  the 
appropriate  junction  temperature  to  update  the  transport  and  detrapping  parameters. 

In  this  paper,  a physics-based  model  is  reported  to  determine  the  power  performance  of 
GaN-MESFETs  by  considering  dispersion  in  the  output  resistance  and  transconductance  due 
to  traps  while  accounting  for  the  thermal  effects.  The  present  treatment  will  allow  the 
optimization  of  GaN-based  MESFETs  used  in  high  power  and  high  frequency  applications. 

Analysis 

In  the  present  analysis,  the  intrinsic  and  the  trap  related  circuit  parameters,  as  shown  in 
Fig.  1,  are  obtained  from  a physics  based  analysis.  Temperature  and  field  dependent  transport 
parameters  are  obtained  from  an  ensemble  Monte  Carlo  simulation  [8,  9].  Following  the 
treatment  reported  by  Golio  et  al.  [12]  the  effects  of  traps  are  incorporated  through  the 
parameters  Css,  g"m  and  R^.  These  parameters  are  obtained  once  the  underlying  physical 

processes  governing  trap  dynamics  are  formulated.  Intrinsic  circuit  parameters:  Cgs,  Cgdi  Sm  ’ 
Ri  and  are  obtained  from  conventional  small-signal  MESFET  analysis  in  the  absence  of 
traps  and  considering  the  effects  of  velocity  saturation.  The  other  circuit  parameters  are 
obtained  from  the  reported  experimental  data  and  are  as  follows:  Rg  = 6£l,Lg  = 0.055  nH,  Rd 
= 70  Q,  Ld  = 0.307  nH,  Rs  = 90  Ls  = 0.027  nH  and  C*  = 0.040  pF  [1 1, 12]. 

A thermal  simulator  is  incorporated  to  calculate  the  operating  temperature  for  a given 
power  level  by  solving  Laplace’s  equation.  The  analysis  proceeds  with  the  initial  guess  of 
junction  temperature  to  estimate  mobility,  carrier  velocity  and  critical  electric  field  for  a 
given  channel  electric  field  to  determine  the  drain  current  for  given  drain  and  gate  bias. 
Estimating  the  dissipated  power  which  is  related  to  the  product  of  drain  voltage  and  current 
Laplace’s  equation  is  solved  to  obtain  the  updated  temperature.  Transport  parameters  and 
critical  electric  field  are  updated  to  determine  the  drain  current  and  the  process  continues  till 
a consistent  solution  in  terms  of  drain  current  and  temperature  is  obtained.  The  temperature 
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dependent  transconductance  and  output  resistance  dispersion  frequencies  are  obtained  once 
self-consistency  is  achieved. 

At  any  given  frequency  the  overall  intrinsic  output  conductance,  g^,  and 
transconductance,  gm,  are  obtained  by  analyzing  the  output  subcircuit  of  Fig.l.  For  large- 
signal  analysis,  gm,  gds  and  Cgs  are  considered  nonlinear  functions  of  vg  while  Cgd  is  nonlinear 
function  of  vdg.  The  nonlinear  functions  are  approximated  up  to  second  order  term,  p = po  + 
Pi  v + p2  v , where  p represents  gm,  gds,  Cgs  or  Cgd  and  v represents  vg  or  vdg. 

- Vdg  + 


LS  Rg  csd  Rd  Ld 


Fig.l.  GaN  MESFET  model. 

Volterra  series  technique  is  used  to  compute  output  power  and  nonlinearities  at  RF.  For 
Vol terra  series  analysis,  the  circuit  shown  in  Fig.l  is  terminated  in  source  impedance  Zs  and 
load  impedance  ZL  across  the  gate-source  and  drain-source  terminals,  respectively.  Defining 
ports  1 through  5 across  nonlinear  elements  Cgs,  gm  and  g*,  Cgd,  load  and  source  respectively, 


the  elements  of  5 x 5 system  matrix  Y are  expressed  as  follows: 

YM=yRi+jocgs 0,  a) 

yi,2  W = W = Yxi  W = -i /Rf , (2) 

Y2,i(a>)=-yRi+gmo , (3) 

^1.4  W = YA,l  W = ^1,5  W = Yu  {co)  = 0 , (4) 


^2.2 (^)-  l/*f  + [{V(^4  + j&Ld)+  \/(zs(co)+  Rg  + jcoLg  )}  1 + Rs  + jcoLs}'  + j(!)Cds  + gds0 , (5) 
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^.4  (^)= [Mz*(®)+ Rg + j^yy^s + +r„+ ja)LdY + yzL{a>), 

Y5,s  M = l{l/fe  + joh ) + Vfo  + j<»Ld 

+ Zs(a>)  + Rg  + j coLg  ]"* 

*2,5  (®)  = *5,2  (®)  = “*5,5  M 

[(Rd  + jcoLd)/(Rs  + j(dLs  +Rd+  jcoLd )] 

y2A  M« Yu  W=lvz»  - r4»J 

/?s  + jmL,  + zs(<y)+i?g  + jgjL?  _ 
y3» = r4» = [y4>)-i/z>)J 
h±j®k » 

Rs  + jcoLs  + Zs(a))+Rg  + jo)Lg 

n,3W=n.2W=v^-^W* 

*3,3  (®)  = *5,5  (®)  V*,  + j®Cdg 0 ’ 

*3,5  M = *5,3  W = *5,5  M . and 


(6) 

(7) 

(8) 
(9) 


(10) 

(ID 

(12) 

(13) 

(14) 


With  an  applied  signal  of  amplitude  Vs\  at  frequency  ah,  voltages  across  ports  1 through  4 
are  determined  using  [14]: 


[vP  («>,)]=  -k.,., 


(15) 


where  p = 1,...,4.  The  second-order  voltages  appear  across  ports  at  mixing  frequencies,  cih,k  - 
(Q\  + Oh  when  two  tones  are  applied  with  amplitudes  Vsi  and  VS2  at  frequencies  (o i and  Oh, 
respectively.  Second-order  port  voltages  are  calculated  by  applying  nonlinear  currents 
through  each  nonlinear  port.  The  nonlinear  currents  are  evaluated  using  first-order  port 
voltages  due  to  individual  tones  and  p\  coefficients  of  nonlinear  elements  [14].  With  Y 
matrix  evaluated  at  mixing  frequencies,  the  second-order  port  voltages  are  given  by, 


K (®i*®a)]=  4W>,-  J"'  Luf- P = h-4 . 


(16) 


Similarly,  third-order  port  voltages  due  to  tone  amplitudes  Vs\,  Vsi  and  V*3  and 
frequencies  ah,  Oh  and  ah,  are  calculated  by  applying  current  sources  due  to  nonlinear 
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coefficients  p\  and  P2  and  first  and  second-order  port  voltages.  The  fundamental  component 
of  the  output  power  is  given  by, 

K.  =0.5|Fi(I»1)/ZI(©,)|!Re[Zl(ffl1)].  (17) 


Results  and  Discussion 

Fig.2  shows  calculated  DC  and  pulsed  I-V  characteristics  incorporating  thermal  and 
trapping  effects  for  a 0.3/tfn  x 100/an  GaN  MESFET  [10],  Experimental  results  are  plotted 
on  the  same  figure  to  show  good  agreement  [10].  The  MESFET  structure  was  grown  on 
sapphire  with  a 2000  A n-GaN  active  layer  with  doping  concentration  of  2.7  x 1017  cm'3  [10]. 
In  the  presence  of  light,  the  DC  I-V  measurements  do  not  show  any  current  collapse  as  the 
electrons  captured  by  traps  located  at  the  channel-buffer  interface  and  the  surface  (between 
gate-drain  region)  have  sufficient  time  to  be  released.  With  pulsed  input  signals  electrons 
captured  by  the  buffer  traps  form  a depletion  region  in  the  channel  at  the  channel-buffer 
interface.  Besides,  surface  traps  also  capture  electrons  and  form  a virtual  gate  [15]  which 
causes  drain  current  to  decrease  for  a given  drain  bias.  The  effect  of  surface  traps  is 
incorporated  by  considering  an  additional  negative  gate  potential  due  to  the  trapped  electrons 
in  the  region  between  gate  and  drain. 


Fig.2.  Calculated  (solid  lines)  and  measured  (symbols)  I-V  characteristics  for  0.3//m 

x 100/an  GaN  MESFET  [10]. 
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Input  power  (dBm) 

Fig.3.  Variation  of  output  power  as  a function  of  input  power  for  Lxl00//m  GaN 

MESFET  at  4GHz. 


Fig.3  shows  the  variation  of  the  output  power  as  a function  of  input  power  for  0.3//m  and 
0.75//m  gate  length  devices  operating  at  4GHz.  The  bias  voltages  are  Vgs  - -2V  and  Vos  = 
30V.  Using  isothermal  (300K)  calculations  maximum  output  powers  of  27dBm  and  19dBm 
are  obtained  with  0.3/rni  and  0.75//m  gate  length  devices,  respectively.  For  the  same  gate 
lengths  by  taking  into  account  thermal  effects  the  maximum  output  powers  decrease  to 
22dBm  and  17dBm,  respectively.  Higher  output  power  is  obtained  for  a given  input  power  in 
shorter  gate  length  devices  due  to  higher  transconductance  resulting  from  higher  differential 
mobility.  The  temperature  and  electric  field  dependent  differential  mobility,  defined  as  the 
ratio  of  velocity  to  electric  field,  for  the  0.3/nn  gate  length  device  is  given  by,  jUn(E,T)  = 
(1705.2  - 30.912s  + 0.1564E2)  + (-3.35  + 0.075E-  4xlO^F2)T+  (0.0015  - 4xl0_5£:  + 2x10“ 
7Zs2)72,  where  temperature  T is  in  K and  channel  electric  field  E is  in  kV/cm.  Mobility 
expression  for  the  0.75//m  gate  length  device  is,  /Jn(E,T)  = (1914.4  - 36.18E  + 0.18552s2)  + (- 
4.30  + 0.092T  - 0.000532s2)r  + (0.0032  - 8xl0_5£  + 5xl(T7E1)fl.  The  determination  of  the 
temperature,  electric  field  and  size  dependence  of  mobility  proceeds  by  carrying  out  an 
ensemble  Monte  Carlo  simulation  by  taking  into  account  scatterings  due  to  polar  optical 
phonon,  acoustic  phonon,  equivalent  and  non-equivalent  intervalley  scattering,  impact 
ionization,  ionized  impurity,  alloy,  interface  and  self-scattering  [8,  9].  For  a given  power 
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dissipation,  shorter  gate  length  devices  operate  at  higher  temperatures  and  their  mobility  and 
saturation  velocity  are  greatly  reduced  from  their  corresponding  room  temperature  values. 
Though,  at  a given  temperature  and  electric  field  shorter  gate  length  devices  have  a higher 
differential  mobility  as  compared  to  devices  with  longer  gate  lengths.  The  higher  temperature 
in  shorter  structures  is  due  to  higher  dissipation  power  density.  As  a result,  output  power 
decreases  significantly  due  to  device  heating  in  shorter  structures. 


Fig.4.  Variation  of  power  gain  as  a function  of  input  power  for  LxlOO/an  GaN  MESFET 

at  4GHz. 


In  Fig.  4,  a similar  result  showing  gain  compression  due  to  device  heating  is  depicted. 
The  gain  compression  due  to  thermal  effects  for  the  0.75/an  gate  length  device  is  0.75dB. 
Due  to  higher  operating  temperature  the  gain  compression  for  the  0.3/an  gate  length  device  is 
2.2dB.  Similar  effects  in  power  added  efficiency  (PAE)  is  shown  in  Fig.5.  Shorter  gate  length 
devices  have  better  PAE  due  to  higher  power  gain.  PAE  compression  due  to  thermal  effects 
is  less  in  longer  gate  length  devices  and  follows  the  trend  observed  in  output  power  and  gain 
compression. 
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0 5 10  15  20 

Input  power  (dBm) 

Fig.5.  Variation  of  power  added  efficiency  (PAE)  as  a function  of  input  power  for 

LxlOO/zm  GaN  MESFET  at  4GHz. 


Conclusions 

RF  power  performances  of  GaN  MESFETs  are  analyzed  considering  trapping  and 
thermal  effects.  The  DC  to  RF  dispersions  of  transconductance  and  output  resistance  due  to 
detrapping  effects  and  transport  parameter  variations  with  channel  temperature  are 
incorporated.  Volterra  series  technique  is  used  to  analyze  device  power  performance  and 
nonlinearities.  Short  channel  devices  demonstrate  higher  output  power,  power  gain  and  PAE, 
however,  better  thermal  stability  is  obtained  in  long  channel  devices. 
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